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The interplay between structure, magnetism and superconductivity in single crystal 
Ba(Fei_ ;c Co :! ;)2As2 (x=0.047) has been studied using high-resolution X-ray diffraction by monitor- 
ing charge Bragg reflections in each twin domain separately. The emergence of the superconducting 
state is correlated with the suppression of the orthorhombic distortion around T c , exhibiting com- 
petition between orthorhombicity and superconductivity. Above Ts, the in- plane charge correlation 
length increases with the decrease of temperature, possibly induced by nematic fluctuations in the 
paramagnetic tetragonal phase. Upon cooling, anomalies in the in-plane charge correlation lengths 
along a (£ a ) and 6 axes (£;,) are observed at T s and also at T B indicative of strong magnetoelastic 
coupling. The in-plane charge correlation lengths are found to exhibit anisotropic behavior along and 
perpendicular to the in-plane component of stripe-type AFM wave vector (101)o below around Tn. 
The temperature dependence of the out-of-plane charge correlation length shows a single anomaly 
at Tn, reflecting the connection between Fe-As distance and Fe local moment. The origin of the 
anisotropic in-plane charge correlation lengths £ a and £(, is discussed on the basis of the antiphase 
magnetic domains and their dynamic fluctuations. 

PACS numbers: 74.25.Ha, 74.70.Xa, 75.30.Fv, 75.50.Ee 



INTRODUCTION 

In the recently discovered iron-based superconductors 
PJ [2], the superconducting temperatures are found to 
be in close proximity to an antiferromagnetic (AFM) 
and a tetragonal-orthorhombic (T-O) structural transi- 
tions. It turns out that the suppression of both the 
AFM and the T-0 structural transitions, by doping 
or by pressure, eventually induces superconductivity. 
These phenomena beg the question about the role of 
spin and lattice degrees-of-freedom in the emergence of 
superconductivity [2J [3]. Several theoretical descriptions 
have been proposed to interpret the relationship between 
the structural and AFM transitions on the basis of or- 
bital ordering or by introducing an intermediate spin- 
nematic phase resulting from an effective — J 2 local- 
spin model or from an itinerant model both with equiv- 
alent consequences. [3HS] Although different in details, 
these descriptions emphasize the importance of the mag- 
netoelastic coupling in driving the two transitions simul- 
taneously or separately. Canoei al. [5] studied the inter- 
play between the elastic and spin degrees-of-freedom in 
iron pnictide superconductors using a Ginzburg-Landau 
approach, indicating that the magnetoelastic coupling 
can bring about the particular features of the structural 
and magnetic transitions in these systems including the 
emergence of the collinear stripe-type AFM ordering. Re- 
cently, a microscopic study [7] of a simple symmetry- 
allowed model Hamiltonian demonstrated that due to 
the effect of magnetoelastic coupling, the considerable 
orthorhombic clastic softening is caused by critical spin 
fluctuations present in the system before magnetic or- 



der occurs. This may explain why the AFM transition 
is often preceded by the T-0 structural transition. It 
should be pointed out that this picture is similar to 
the nematic phase model[31|S]. To date, there are very 
few experimental reports on the magnetoelastic effect 
in iron pnictides and such reports investigated the role 
that magnetoelastic effect plays in the structural and 
magnetic transitions under the application of external 
driving forces. Magnetoelastic effects have been demon- 
strated by applying pressure to CaFe2As2 and inducing 
an O-T and AFM-to-nonmagnetic transitions [5], or by 
applying shear stress to BaFe2(Asi_ a; P; C )2 that shifts the 
magnetic transition and superconducting critical temper- 
atures significantly. [9] Therefore, it is of interest to study 
the possibly intrinsic magnetoelastic effect in iron-based 
superconductors, without introducing any external driv- 
ing force. 

The Co-doped BaFe 2 As2 system exhibits a rich phase 
diagram with a complex interplay between the struc- 
tural, magnetic, and superconducting phases. [I0HT3] In 
the parent BaFe2As2 compound, the AFM ordering tran- 
sition at T N coincides with a T-0 structural transition 
at T s . Upon doping both transitions gradually separate, 
such that T s > T N , accompanied with the appearance of 
superconductivity above x = 0.03 in Ba(Fei_ a; Co 2 ;)2As2. 
It is interesting to point out that the orthorhombic dis- 
tortion 6 in Ba(Fei_ :c Co : j : )2As2 with lower Co content 
(x = 0.018) shows one clear anomaly at the magnetic 
transition temperature T N , but it is absent at Tn for 
x = 0.047 superconductor with intermediate Co con- 
tent. For x higher than ~0.066, both the magnetic 
and structural transitions are completely suppressed and 
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FIG. 1: Schematic illustration in reciprocal space of (a) un- 
twinned crystal in the (HK0)t plane at the high-temperature 
tetragonal phase and (b) at the orthorhombic phase with 
uniformly rotated and separated anisotropic twin-domains in 
(hkQ)o plane. The arrows show typical scans performed at 
Bragg reflections. 

superconducting transition is the only transition ob- 
served. Only in an intermediate composition region of 
0-03 < x < 0.066, does Ba(Fei_ :E Co K )2As2 exhibit a 
coexistence of superconductivity, O-structure and AFM 
phases providing potential candidates to investigate the 
effects of the magnetoelastic coupling. However, the ten- 
dency of these crystals to form twinned orthorhombic 
domains has hampered definitive determination of inher- 
ent features of the intermediate phase between T s and 
T N where the presumed nematic phase exists. Therefore, 
there have been extensive efforts to de-twin these crys- 
tals [14] to establish the underlying electronic, structural, 
and magnetic anisotropics that characterize this inter- 
mediate phase. Motivated by these issues, we set out 
to investigate anisotropic features of the crystal struc- 
ture and domain formation over a wide range of tem- 
peratures using high resolution x-ray diffraction meth- 
ods that reveals the intrinsic magnetoelastic coupling in 
the Ba(Fei_ x Coa;)2As2 (x — 0.047) superconductor. The 
high resolution allows us to separately monitor Bragg 
reflections of different orthorhombic twin domains and 
study their temperature evolution, as has been done re- 
cently on CeFeAsO[J5]. 

EXPERIMENTAL DETAILS 

The Ba(Fei_ :E Co :E )2As2 (x = 0.047) crystal was 
grown using a self-flux solution method as described 
previously [HI [TB]. The crystal has been character- 
ized by X-ray diffraction, neutron scattering, magneti- 
zation, resistivity and heat capacity, identifying reported 
three transitions, structural at T s ~ 60 K, magnetic at 
T N ~ 47 K and SC at T c ~ 17 K.pU US] A plate- 
like piece of crystal with its c-axis perpendicular to its 
surface was chosen for investigations by high-resolution 
x-ray scattering technique using the six-circle diffrac- 



tometer of the 6-ID-B beamline at the Advanced Pho- 
ton Source (APS), Argonne National Laboratory (X-ray 
energy kept at 8 keV) . The scattering geometry of our ex- 
periment is similar to that used previously by Li et al. [15] . 
We use orthorhombic indices {hkl)o at all temperatures 
so that the tetragonal (HKL)i indices are provided in 
terms of the twin domains in the orthorhombic structure 
with the following conversion (H + K,H — K, L)o and 
(H — K, H + K, I/)o- The crystal was mounted at the end 
of the cold-finger of a Displex cryogenic refrigerator with 
access to (00Z)o and high index (hkt)o Bragg reflections. 
Flux intensity on the sample was optimized to eliminate 
beam heating effects of the sample while maintaining a 
reasonable signal to noise ratio. To achieve that, slit 
setups and attenuations yielding Bragg reflection inten- 
sities that scaled with beam attenuations were chosen. In 
this regard, it should be noted that the low thermal con- 
ductivity in the SC state required a significant increase 
of beam attenuation to prevent sample-heating during 
measurements. These considerations limited the choice 
of setups, i.e., resolution, but as discussed below, by ade- 
quate analytical tools we captured the intrinsic behavior 
(i.e., charge correlation lengths, in particular) of this sys- 
tem. As demonstrated on the CeFeAsO|15j. the twin do- 
mains are uniformly rotated and separated in reciprocal 
space by a microscopic shear-angle enabling the charac- 
terization of each domain. Fig. 1 shows schematically a 
limited in-plane reciprocal zone with the (110)t of the 
untwinned crystal (a) that transforms to the (200) o and 
(020)o in the orthorhombic symmetry notation (b). The 
misfit angle between the two domains is determined by 
the orthorhombic distortion ( 8 = ) . [TS] The arrows in 
Fig. [T] show typical scans performed at Bragg reflections. 
In the present study, we monitored the (208)o/(028)o 
systematically (corresponding to (118)t at high temper- 
atures) and also the (008)o, as a function of temperature. 

RESULTS AND DISCUSSION 

Fig. 2(a) shows the evolution of the orthorhombic- 
ity as the temperature decreases with a sharp splitting 
of the (118) T into the (208) o and (028) o at T s = 60 K. 
With further decrease in temperature, the orthorhombic- 
ity increases without displaying an anomaly at T" N = 47 
K (see the inset of Fig. 2(a)), which is consistent with 
the report by Kim et al. [THJ. A slight decrease in the 
splitting is observed below the SC temperature T c , indi- 
cating that superconductivity and the orthorhombic dis- 
tortion are coupled. [10 . Note that for higher Co substi- 
tution, the suppression of orthorhombic order parame- 
ters becomes larger. For example, the orthorhombic dis- 
tortion in x = 0.063 is completely suppressed and the 
reentrant transition to tetragonal structure occurs below 
T c . When x increases to 0.066, the orthorhombic dis- 
tortion vanishes and no T-0 transition is observed. [10] 
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FIG. 2: (color online) The temperature dependence of peak 
positions extracted from (a) (/i08)o (h domain, circles) and 
(0fc8)o (k domain, squares) scans, and (b) the I scan for (00£)o 
reflection. The inset of (a) shows the orthorhombic distor- 
tion S as a function of temperature. The inset of (b) shows 
the zoomed view on I scan for (20£)o reflection. Linear and 
quadritic fit to the data above Ts are included as solid lines. 
The dashed lines in (a) and (b) indicate the locations of struc- 
tural, magnetic and superconducting transition temperatures 
Ts, In, and T c . 



From Fig. 2 (a), we point out that the in-plane lat- 
tice constant shrinks linearly in the tetragonal phase of 
Ba(Fei_ x Co a; )2As2 with decreasing temperature at a rate 
of 1.4 x 10~ 5 A/K per unit cell (linear thermal expansion 
parameter a ~ 2.5 x 10~ 6 /K). By contrast, the thermal 
expansion along the c-axis is weakly quadratic (Fig. 2(b)) 
with no abrupt anomaly at T s and displays a deviation 
from the quadratic form near T N (a weak minimum is 
observed in c-axis lattice parameter for (20Z)o reflection 
as shown in the inset of Fig. 2(b)). The effect of mag- 
netic transition at T N on the lattice parameter c implies 
a coupling between them, which will be discussed below. 

Fig. 3(a) shows representative in-plane scans along h 
for the (208)o Bragg reflection at 150 K (above T s ) and 
46 K (below T s ). It is clear that the full width at half 
maximum (FWHM) of the Bragg peak is broader be- 
low T s . Variation in the line broadening can be due to 
changes in coherence length, domain size, mosaic distri- 
butions, and more likely a combination of all three. To 
obtain quantitative evaluation of peak line-widths, these 
and other scans were initially modeled as a Gaussian, 
a Lorentzian, their linear combination, or Pseudo-Voigt 



line-shapes, but none of these lineshapes yielded satisfac- 
tory agreement with the data. We therefore adopted a 
standard convolution method by systematically folding a 
Gaussian resolution function and a Lorentzian function 
that reflects an exponentially decaying charge (chemical) 
coherence length as follows, 



G(q')L(q - qi)d q > 



(1) 



where G(x) 



-{\n2)(x) 2 /ui 2 



and L(x) 



ujV tt In 2 

= , , c '-, , 2 , such that 2uj and 2v are the FWHM of the 
Gaussian and Lorentzian functions in reciprocal space, 
respectively. Our resolution was high enough to resolve 
the twin domains separately, i.e., optimizing peak inten- 
sity of the (208) o from one twin domain and that of the 
(028)o of the other domain required sample rotation be- 
tween the two peaks as has been done for CeFeAsO|15j. 
While the FWHM of the Lorentzian function (2v) repre- 
sents the intrinsic width n of the sample, we should note 
that the Gaussian function in Eq. has two contribu- 
tions: one from geometrical setup (i.e., incident and scat- 
tered beam divergence) and the other from the mosaic 
spread of the studied crystal (see a detailed discussion in 
resolution function in Ref. [II]). It is by now well estab- 
lished that the mosaic spread of typical pnictides under- 
going shear induced O-T transition exhibit mosaic spread 
changes due to stresses during the transitions. We there- 
fore attribute the temperature dependence of the Gaus- 
sian width as arising primarily from the variation in the 
mosaic distributions. Fig. 3(b) shows the FWHM of the 
Gaussian function for the (208)o along /iasa function of 
temperature indicating anomalies that can be related to 
the stresses introduced by the structural and magnetic 
transitions in the system. 

The temperature evolutions of the intrinsic width 
k (= 2v) extracted from the Lorentzian functions of 
(208)o and (028)o Bragg peaks are shown in Fig. 4(a). 
It is evident that anomalies are observed at both T s 
and T N . Similar observation by using high-resolution 
x-ray diffraction was reported in different systems 
TbVi_ a; As :I ;04 (x=0 and 1) [IS], where the intrinsic 
width of the Bragg peaks shows a clear peak at their 
T c that is reminiscent of A anomalies in the heat capac- 
ity. Note that the intrinsic width k in unit of -A -1 , i.e., 
the FWHM of the Lorentzian profiles from X-ray scatter- 
ing in reciprocal space corresponds to the inverse charge 
correlation length £ [T51 ITBH213] : 



(2) 



Thus, the temperature dependence of the in-plane 
charge correlation lengths along and normal to the 
in-plane component of AFM propagation wavevector 
(101)o, i-e., longitudinal charge correlation length £ a 
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FIG. 3: (color online) (a) A representative longitudinal scan of 
Q=(/i08)o at 150 K (above T s )) and 46 K (below T s )). The 
Q values and intensity were normalized for comparison. It 
can be seen that the peak widths broaden significantly below 
T s . (b) The temperature evolution of FWHM of the Gaussian 
function (2ui) for the longitudinal scan of h domain. 
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FIG. 4: (color online) The temperature dependence of (a) 
the intrinsic width « of the longitudinal scans of h and k 
domains obtained from the Lorentzian function and (b) the 
in-plane charge correlation lengths along a axis (£ n ) and b 
axis (£(,). The inset of (a) shows a detailed view of k around 
Tn and T s - The anisotropy of the correlation length rj^ as 
a function of temperature is shown in the inset of (b). The 
dashed lines mark the locations of structural and magnetic 
transition temperatures Ts and T s . 



(along the AFM bond direction) and transverse charge 
correlation length ^ (along the ferromagnetic bond direc- 
tion), can be derived from the intrinsic widths of (208) o 
and (028)o Bragg peaks. As illustrated in Fig. 4 (b), 
both £ a and £b show two clear peaks at T s and T N , 
respectively. With the decrease of temperature to T s , 
the in-plane charge correlation lengths increase gradu- 
ally, followed by a rapid decrease below T s . When the 
temperature approaches T N , the in-plane charge corre- 
lation lengths increase again. It is worth emphasizing 
that the Bragg reflections used to extract these data 
are not allowed by the symmetry of the magnetic struc- 
ture of the ordered iron moments and they are strictly 
the result of charge (nuclei) ordering. Interestingly, 
the intrinsic width k and charge correlation length ob- 
tained from the high-resolution X-ray data show a clear 



anomaly at the magnetic transition temperature, sug- 
gesting these charge Bragg reflections are sensitive to the 
spin-structure and fluctuations. This is presumably due 
to the strong magneto-elastic coupling that exerts sec- 
ondary effects on charge correlations, domain formation 
and their shape. 

Below around T N , we notice that the charge correlation 
lengths along and normal to the in-plane component of 
AFM propagation wavevector (101)o are different (£ a > 
£b), displaying an anisotropic behavior. The anisotropy 
in charge correlation length is defined as [21] 



f 2 - f 1 

So ?b 



(3) 



rjc=0 indicates isotropic correlations (£ a = £&)i whereas 
i]^=l (ry^=-l) corresponds to the extreme case of £ a ^ 
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(£o £b) for structural domains consisting of long lin- 
ear stripes. Based on this equation, we have derived the 
anisotropy of the charge correlation length as a function 
of temperature, as shown in the inset of Fig. 4. The 
anisotropy is most pronounced below T N , with values 
ranging from 0.1 to 0.8 as the temperature is lowered. 

A possible scenario to interpret the anisotropic charge 
correlation lengths along a and b below around T$ in 
x = 0.047 may be related to the presence of the antiphase 
magnetic domains 22 - 24 . Mazin and Johannes [22] first 
proposed that antiphase domains and their dynamic fluc- 
tuations are central for understanding the high- T c fer- 
ropnictides. Very recently, Li et.al [21] observed surface- 
pinnned antiphase domains in BaFe2As2 using high- 
resolution scanning tunneling microscopy. Since the en- 
ergy differences between the AFM stripe magnetic struc- 
ture and other AFM patterns are small [22], it is highly 
possible that many antiphase magnetic boundaries are 
formed. The antiphase domains are pinned at T < T N , 
and show dynamic fluctuations in the region of T N < T 
< T s . There are two kinds of simple antiphase domains 
(labeled A and B) with boundaries along a and b axes, as 
shown in Fig. 5. Due to the same magnetoelastic inter- 
actions that lead to a difference in the ferromagnetic and 
AFM bond lengths in the orthorhombic structure, we 
propose that the formation of such antiphase domains 
are accompanied by elastic distortions at their bound- 
aries. The antiphase boundaries along a axis influence 
the magnitude of transverse charge correlation length 
whereas the antiphase boundaries along b axis affect the 
longitudinal charge correlation length £ Q . Differences in 
the density of antiphase boundaries in the two directions 
eventually leads to the anisotropy in £ a and below Tn- 
The dynamic fluctuations of the antiphase domains may 
be responsible for the anisotropy in charge correlation 
lengths in the small temperature region above T N (but 
lower than Ts), as shown in the inset of Fig. 4 (b). 

Based on such antiphase magnetic domain scenario and 
the Ji — J 2 model, the anisotropic charge correlation 
lengths can be used to estimate the ratio of the mag- 
netic exchange parameters J\ and J 2 . As shown in Fig. 
5, the energies per-spin (S) for forming these two kinds 
of antiphase boundaries are given by 



E A = (2J 2 ^J 1 )S 2 



E B = (2J 2 + J^S 2 



(4) 



(5) 



Since J\ is AFM, the number of antiphase domain 
boundaries Na arid Nb should scale inversely with the 
magnetic energy of the domain wall so that 
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FIG. 5: (color online) Schematic pictures of two kinds of 
the antiphase magnetic domains, with boundaries along (a) 
a axis and (b) 6 axis, respectively. The antiphase domains 
are pinned at T < T H , but show the dynamic behavior in the 
region of Ts < T < T s - The pink lines show the antiphase 
boundaries. J ± represents the nearest-neighbor exchange cou- 
plings along a or b directions, whereas J 2 represents the next- 
nearest-neighbor exchange couplings. 



The charge correlation length scales inversely propor- 



tional to the number of boundaries, i.e., £ a oc 

& = Na = 2 J 2 + Ji 
& Nb 2J 2 - J x 



N L 



Thus, 



(7) 



(6) 



In Fig. 4 (b), we observe £ a /6> w 3 at low temper- 
atures and from Eq. (6) we can get J 2 J\, which 
is consistent with the previous calculations or experi- 
ments on other iron pnictides, such as LaFeAsO [5S] and 
Ba(Fei_ x Co :r )2As2 (x = 0.074)[2H|. This is also reason- 
able for producing a stripe-type AFM structure that re- 
quires J 2 > 

It is worthwhile noting that the in-plane charge cor- 
relation lengths show gradual changes at temperatures 
significantly above T s , probably due to magnetic fluc- 
tuations that are known to persist above T s in similar 
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FIG. 6: (color online) Temperature dependence of (a) the in- 
trinsic width ft of the / scan for (008) o obtained from the 
Lorentzian function and (b) the out-of-plane charge correla- 
tion length along c axis (£ C )-The dashed line indicates the 
locations of T« and TV 



pnictides[271 [3H] - We point out that this feature may 
support the nematic model. In the nematic model[3l[6], 
the nematic order coincides with the structural tran- 
sition with the notion that the driving force for the 
T — O transition is not elastic in origin but magnetically 
driven by Ising-like interpenetrating AFM domains [3J [7] . 
Nematic (magnetic) fluctuations remain at higher tem- 
perature above T s , which has been suggested by var- 
ious techniques, such as susceptibility anisotropy [25] . 
shear modulus [3j, inelastic neutron scattering! 30j, and 
anisotropic in-plane resistivity [H]. 

We now turn to a discussion of the c-axis charge cor- 
relation length. Fig. 6(a) and (b) show the temperature 
dependence of the intrinsic width of the (008) o Bragg 
reflection and the corresponding out-of-plane charge cor- 
relation length along c axis (£ c ), respectively. A single 
anomaly at the magnetic transition temperature T N and 
a gradual change above and through Ts, are observed. 
The absence of sharp anomaly in £ c suggests that atomic 
distortions resulting from the T-0 structural transition 
mainly occur in the ab plane. This is consistent with the 



fact that the in-plane lattice parameters a and b change 
significantly, but c changes weakly around T s , as can also 
be seen from Fig. [2j Both the out-of-plane charge corre- 
lation length £ c and lattice parameter c show an anomaly 
at T N , showing close correlation between the AFM mag- 
netic transition and the modification of structure along c 
axis. Recent experiments and calculations reveal that the 
Fe local moment is very sensitive to the Fe-As distance 
in iron pnictides. Yin et al. |31j performed density func- 
tional theory (DFT) calculations within the generalized 
gradient approximation and found the Fe-Fe transverse 
exchange coupling is strongly dependent on both the 
AFM symmetry and the Fe-As distance. Belashchenko 
et al. |32j demonstrated that in layered iron-pnictide com- 
pounds, as the Fe-As distance is decreased, the degree of 
itinerancy of Fe moments increase. Moreover, the cou- 
pling between the local moment and the Fe-As distance is 
controlled by strong covalent Fe-As bonding. Recently, 
neutron diffraction studies of CeFeAsi-^P^O [33] and 
DFT calculations [34 demonstrated that a decrease in 
Fe-As distance induces strong hybridization between Fe 
3d and As 4p orbitals, leading to quenched Fe magnetic 
moments. Therefore, the AFM transition at T N is cou- 
pled to the change of Fe-As distance, which leads to the 
anomalies in the out-of-plane charge correlation length 
and lattice parameter c. 

In summary, high resolution X-ray diffraction stud- 
ies on structural Bragg reflections of the SC and AFM 
Ba(Fe 1 _ a; Co :E )2As2 (x=0.047) single crystal reveal sec- 
ondary effect stemming from the magnetic properties 
of the system, which is understood to result from in- 
trinsic and strong magnetoelastic coupling. In addition 
to showing anomalies around the structural and mag- 
netic transitions, the in-plane charge charge correlation 
lengths along a and b axes show anisotropy below around 
T N , which probably results from the effect of antiphase 
boundaries formed along a and b axes. Employing our 
anisotropic charge correlation lengths, we are able to es- 
timate the ratio of J<xj J\ to be around 1 on the basis 
of such antiphase magnetic domain scenario and J\ — 3^ 
model. The out-of-plane charge correlation length £ c and 
lattice parameter c exhibit a single anomaly at T" N , which 
can be associated with the modification of Fe-As distance 
when the AFM transition occurs. Our results also show 
gradual evolution of the Bragg peak widths and in-plane 
charge correlation length above T s , which is presumably 
induced by the nematic magetic fluctuations up to almost 
200 K. 
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